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Dynamic  operation  is  known  as  one  of  the  factors  for  accelerating  chemical  degradation  of  the  polymer 
electrolyte  membrane  in  a  polymer  electrolyte  membrane  fuel  cell  (PEMFC).  However,  little  effort  has 
been  made  dealing  with  the  quantification  of  the  degradation  process.  In  this  investigation,  cyclic  current 
operation  is  carried  out  on  a  fuel  cell  system,  and  the  frequency  effect  of  cyclic  operation  on  chemical 
degradation  is  investigated.  The  dynamic  behavior  of  a  fuel  cell  system  is  analyzed  first  with  the  modi¬ 
fied  Randles  model,  where  the  charge  double  layer  is  modeled  by  three  components;  a  charge  transfer 
resistance  (Rct),  and  two  RC  cells  for  the  Warburg  impedance.  After  calculating  each  parameter  value 
through  exponential  curve  fitting,  the  dynamic  behaviors  of  the  three  components  are  simulated  using 
MATLAB  Simulink®.  Fluoride  release  as  a  function  of  the  frequency  of  cyclic  operation  is  evaluated  by 
measuring  the  concentration  of  fluoride  ion  in  effluent  from  a  fuel  cell  exhaust.  The  frequency  effect  on 
chemical  degradation  is  explained  by  comparing  the  simulated  results  and  the  fluoride  release  results. 
Two  possible  reasons  for  the  accelerated  degradation  at  cyclic  operation  are  also  suggested. 
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1.  Introduction 

The  proton  exchange  membrane  fuel  cell  (PEMFC)  has  been 
considered  a  candidate  for  future  vehicular  and  stationary  appli¬ 
cations  and  has  the  advantages  of  high  current  density,  high 
fuel  efficiency,  and  low  temperature  operation,  although  technical 
problems  remain.  One  of  the  main  issues  limiting  the  application  of 
a  PEMFC  is  membrane  degradation  which  is  mainly  caused  by  three 
factors,  mechanical,  thermal  and  chemical  degradation  [1].  If  it  is 
assumed  that  sufficient  temperature  control  of  a  fuel  cell  system 
can  be  realized,  such  that  thermal  degradation  may  be  neglected, 
mechanical  degradation  is  thought  to  be  the  main  cause  of  early 
failure  [1  ],  and  chemical  degradation  usually  results  in  a  significant 
drop  in  the  performance  before  the  failure  [2]. 

A  number  of  previous  research  efforts  have  shown  that  the 
chemical  degradation  is  mainly  attributed  to  chemical  attack  of 
highly  reactive  oxygen  radicals  formed  in  the  fuel  cell  [1-3]. 
Although  the  mechanisms  describing  how  these  radicals  form  and 
attack  the  membrane  are  still  unclear,  it  is  now  generally  con¬ 
sidered  that  the  formation  of  hydrogen  peroxide  (H2O2)  plays  an 
important  role  in  the  chemical  attack.  The  formation  of  hydro¬ 
gen  peroxide  and  reactive  oxygen  radicals  is  a  multi-step  process 
including  oxygen  gas  crossover  and  is  shown  below  in  Fig.  1  [1,3]. 
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Hydrogen  fluoride  (HF)  is  a  direct  product  of  the  radical  attack  on 
the  membrane.  As  such,  the  rate  of  chemical  degradation  can  be 
evaluated  by  measuring  the  fluoride  content  of  the  effluent  from 
the  fuel  cell  exhaust  [4].  As  fluoride  release  continues,  the  mem¬ 
brane  thickness  decreases,  which  is  called  ‘membrane  thinning’  [4], 
and  eventually,  leads  to  the  generation  of  pin-holes  which  further 
degrade  the  system  by  allowing  hydrogen  to  pass  directly  through 
the  membrane  and  “burn”  with  oxygen  on  the  platinum  catalyst  at 
the  cathode  side  of  the  fuel  cell. 

There  are  several  factors  that  may  accelerate  the  chemical 
degradation.  Severe  degradation  due  to  operation  with  low  humid¬ 
ification  and  operation  with  an  open  circuit  voltage,  but  no  net 
current  has  been  reported  and  investigated  by  multiple  researchers 
[5-7].  A  great  deal  of  research  effort  has  been  focused  on  determin¬ 
ing  the  relation  between  those  factors  and  membrane  degradation, 
but  those  relations  have  not  yet  been  fully  clarified.  Liu  and  Case 
reported  that  dynamic  operation  also  accelerated  the  chemical 
degradation  by  showing  that  there  is  about  a  one  order  difference 
in  fluoride  release  between  cyclic  operation  and  constant  operation 
of  a  PEMFC  [8].  However,  there  has  been  little  further  research  on 
the  effect  of  dynamic  operation  on  chemical  degradation. 

The  goal  of  this  research  is  an  examination  of  the  influence  of 
dynamic  operation  on  accelerated  degradation  of  a  fuel  cell  sys¬ 
tem.  The  PEMFC  used  in  this  research  was  a  Ballard  Nexa  1.2  kW 
fuel  cell  system,  and  dynamic  operation  was  realized  using  a  tran¬ 
sistor  switch  controlled  by  an  ON  and  OFF  pulse  signal  as  shown  in 
Fig.  2.  In  this  figure,  RPL  is  a  parasitic  load  resistor  composed  of  the 
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Fig.  1.  Mechanism  of  chemical  degradation  of  a  PEM  membrane. 


fuel  cell  control  system  and  balance  of  plant  (cooling  fan,  compres¬ 
sor,  solenoid  valves),  and  Ri  is  an  exterior  load  resistor.  Different 
operation  conditions  were  modified  by  changing  the  frequency  of 
the  pulse  signal.  For  the  sake  of  this  work,  the  dynamic  behavior 
of  PEMFC  was  first  investigated.  The  Nexa  system  was  modeled 
by  a  simplified  Randles  model  for  which  parameters  were  eval¬ 
uated  using  curve-fits  and  the  dynamic  behavior  of  the  fuel  cell 
system  as  the  pulse  input  was  varied  was  simulated  using  MAT- 
LAB  software.  The  rate  of  chemical  degradation  during  operation 
with  different  frequencies  was  determined  by  measuring  the  con¬ 
centration  of  fluoride  ions  released  in  the  water  drained  from  the 
fuel  cell  system  using  a  fluoride  ion  selective  electrode.  The  fluoride 
is  usually  released  from  both  the  anode  and  the  cathode,  and  the 
release  rate  is  determined  using  the  effluent  collected  from  both 
sides  [2].  Flowever,  the  product  water  of  the  Nexa  system  is  just 
from  the  cathode  side  because  the  system  has  a  dead  end  structure 
where  the  anode  outlet  is  blocked.  Accumulation  of  nitrogen  and 
water  vapor  in  the  hydrogen  stream  gradually  increases  with  time, 
so  the  anode  is  purged  to  flush  out  inert  constituents  periodically 
[9].  It  is  believed  that  the  fluoride  released  from  the  anode  would 
be  exhausted  with  the  anode  gas  when  it  is  purged,  and  there  is  also 
fluoride  exhausted  in  the  form  of  a  gas  from  the  cathode.  Without 
consideration  of  the  amount  of  these  fluoride  gases,  it  is  very  dif¬ 
ficult  to  determine  absolute  amount  of  the  fluoride  released  from 
fuel  cells.  In  this  research,  the  value  of  the  fluoride  emission  rate 
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Fig.  2.  Experimental  setup  for  dynamic  operation  of  the  Ballard  Nexa  fuel  cell  sys¬ 
tem.  Zt  is  a  total  impedance  of  the  Nexa  system. 


measured  using  the  ion  selective  electrode  does  not  mean  the  abso¬ 
lute  value  of  the  total  HF  emission  from  the  fuel  cell  system.  The 
rate  of  chemical  degradation  is  compared  by  comparing  the  flu¬ 
oride  release  rate  of  the  tests  at  different  operating  frequencies. 
Finally,  the  frequency  effects  on  both  simulated  dynamic  behavior 
and  chemical  degradation  were  compared  and  their  close  relation¬ 
ship  was  shown. 

2.  Experimental  setup 

The  schematic  diagram  of  the  Ballard  Nexa  fuel  cell  system  is 
seen  in  Fig.  2.  To  dissipate  the  fuel  cell  power,  a  0.77  £2  resistor 
with  a  1  kW  power  dissipation  rating  was  implemented  by  parallel 
connection  of  a  1  £2  and  three  1 0  £2  resistors  (FSE-1 000,  Huntington 
Electric  Inc.).  The  duty  ratio  of  the  pulse  signal  was  fixed  at  1 5%  and 
50%  in  which  the  15%  duty  cycle  means  15%  externally  loaded  and 
85%  parasitic  load  only,  and  the  pulse  frequency  was  varied  from 
0.0005  Hz  to  1  kHz.  A  transistor  switch  was  constructed  to  perform 
the  step  loading  and  unloading  of  the  system  using  a  high-current 
complementary  transistor  (MJ11032)  and  a  diode.  The  switch  sig¬ 
nal  and  data  acquisition  were  realized  using  a  standard  DAQ  card 
controlled  using  the  MathWorks  software  package  XPC.  Fluoride 
release  rates  during  operation  of  the  fuel  cell  system  were  evaluated 
by  collecting  effluent  from  the  exhaust  and  measuring  the  fluoride 
ion  concentration  using  a  fluoride  ion  selective  electrode  (Hach, 
Model  No.  51928;  a  concentration  range  for  the  measurement  is 
from  0.01  mgL-1  F  to  saturated  solutions,  and  a  slope  change  for 
the  concentration  is  -58 ±3 mV  per  decade  [10]).  Preparation  of 
a  standard  solution  and  calibration  were  performed  by  the  direct 
measurement  method  as  described  in  its  manual  [10].  The  fluoride 
release  rate  was  calculated  by  measuring  the  fluoride  concentration 
and  condensation  rate  of  the  cathode  product  water. 

3.  Results 

3.1.  Modeling  and  simulation 

The  voltage  and  current  changes  during  cyclic  operation  are 
shown  in  Fig.  3.  The  stack  current  cannot  be  completely  turned  off 
because  there  is  always  a  parasitic  load  current  (i’rpl),  so  the  cyclic 
operation  condition  in  this  research  has  no  open  circuit  condition. 
Furthermore,  under  the  assumption  of  well-controlled  humidifica¬ 
tion  of  the  Nexa  system,  the  only  accelerating  degradation  factor  in 
this  experiment  is  assumed  to  be  the  dynamic  operation. 

The  equivalent  circuit  in  Fig.  4  is  used  to  model  the  fuel  cell 
system  as  described  in  a  previous  effort  [11]  in  which  the  charge 
double  layer  is  modeled  by  three  pairs  of  resistors  and  capacitors 
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Fig.  4.  An  equivalent  circuit  of  the  Nexa  system  modeled  by  the  simplified  Randles 
model.  Vocv,  Rm,  Ru  are  open  circuit  voltage,  membrane  resistance,  and  total  load 
resistance  including  RPL  and  Rs  (switch  resistance). 

components:  a  charge  transfer  resistance  ( Rct),  and  two  RC  cells.  By 
assuming  that  the  oxygen  reduction  reaction  rate  in  the  cathode  is 
much  lower  than  the  hydrogen  oxidation  reaction  rate  and  that  the 
charge  double  layer  is  represented  by  a  pure  capacitance  in  parallel 
with  the  charge  transfer  resistance,  the  equivalent  circuit  of  a  fuel 
cell  can  be  simplified  using  just  one  Randles  model  [11].  The  War¬ 
burg  impedance  of  the  Randles  circuit  can  usually  be  approximated 
by  pairs  of  resistances  and  capacitances  (RC  cells)  in  series,  and 
two  RC  cells  are  used  for  the  Warburg  impedance  in  our  research 
[11].  The  first  RC  cell  (Ri  and  Q)  is  used  to  model  the  time  delay 
related  to  water  diffusion  inside  the  membrane,  and  the  second  RC 
cell  (R2  and  C2)  is  used  to  model  gas  diffusion  in  the  gas  diffusion 
layer.  Values  for  each  of  the  parameters  were  calculated  through 
exponential  curve  fitting.  The  fitting  results  using  three  exponential 
curves  are  shown  in  Fig.  5.  The  voltage  transient  occurring  during 
a  load  implementation  is  termed  a  “down  curve”  (switch  ON)  and 
the  voltage  transient  occurring  during  a  load  removal  is  termed  an 
“up  curve”  (switch  OFF).  If  the  modified  Randles  model  in  Fig.  4  is 
employed,  the  down  and  up  curves  can  be  represented  by  three 
exponential  curves.  Therefore,  six  independent  equations  for  eval¬ 
uating  six  parameters  ( Rct ,  Cdi,  Ri ,  C\ ,  R2,  and  C2 )  can  be  set  up  using 
three  exponential  curves  [1 1  ].  In  fitting  the  up  curve,  the  modified 
voltage  curve  for  fitting  is  a  negative  exponential  curve  converging 
to  zero  with  time.  The  fitted  curves,  three  exponential  curves,  and 
their  sum  are  also  negative,  but  they  are  expressed  as  positive  val¬ 
ues  by  adding  the  stabilized  value  of  the  off  state  voltage  as  shown 
in  Fig.  5(b).  The  curve  fitting  and  calculation  are  performed  by  the 
software  package  MATLAB,  and  the  results  show  that  the  parame¬ 
ter  values  calculated  from  the  down  curves  are  not  same  as  those 
of  the  up  curves.  The  noticeable  reason  for  this  difference  can  be 
found  in  compressor  operation  of  the  Nexa  system.  The  fitting  range 
for  the  down  curve  is  performed  for  the  data  collected  prior  to  the 
acceleration  of  the  compressor,  but  the  range  for  the  up  curve  is  at 
the  already  accelerated  compressor  speed  [11]. 

Since  the  down  curve  and  up  curves  are  different  from  each 
other,  the  parameter  values  calculated  from  these  two  curves  are 
also  different.  Each  parameter  value  calculated  from  the  experi¬ 
mental  data  of  50%  duty  ratio  and  0.2  FIz  under  the  assumption  of 
neglecting  the  current  to  the  parasitic  load  (iRPL)  was  summarized 

Table  1 

Parameter  values  of  Randles  model. 


time  /  s 

Fig.  5.  Curve  fitting  using  three  exponential  curves  of  (a)  down  curve  and  (b)  up 
curve. 

in  Table  1.  The  simulation  was  conducted  twice  for  the  down  and 
up  curve  differently  using  MATLAB  Simulink®  as  shown  in  Fig.  6, 
and  the  results  were  plotted  in  Fig.  7. 

Due  to  the  difference  of  parameter  values  between  down  and  up 
curves,  it  is  difficult  to  reproduce  the  continuous  voltage  response 
to  the  high  frequency  pulse  input  especially  in  cases  with  asym¬ 
metric  duty  ratios  such  as  15%.  The  simulation  was  performed  at 
50%  of  duty  ratio  for  checking  the  possibility  to  reproduce  the  volt¬ 
age  profile  at  high  frequency  using  the  parameter  values  calculated 
from  low  frequency  data  of  0.2  Hz.  In  Fig.  8,  it  can  be  seen  that 
the  simulation  results  compare  well  with  experimental  results  of 
various  frequencies.  For  the  purpose  of  the  continuous  plot,  all 
starting  points  of  up  curves  were  brought  into  the  end  points  of 
down  curves.  The  current  to  the  parasitic  load  slightly  changes  in 
terms  of  the  frequency,  so  RLt  should  be  modified  for  simulating 
high  frequency  response.  RLt  at  each  frequency  in  Fig.  8  was  cal¬ 
culated  using  data  between  4  s  and  5  s  of  the  experimental  curve 
where  the  voltage  is  locally  stabilized. 

In  cyclic  operation  of  the  Nexa  system,  a  period  mismatch  is 
shown  between  the  input  signal  and  compressor  operation.  At  low 
frequencies,  the  compressor  speed  is  low  when  the  switch  is  off 


Vocv  (V) RmW_ cdl  (F) Rct  (Q) Ki(fl) Q(F) R2  (Q) C^F) Ru  (Q) 

Down  curve  41.82  0.0908  0.0168  0.0718  0.145  1.668  0.213  0.0631  0.825 

Up  curve  41.82  0.0866  0.0213  0.0807  0.177  1.442  0.187  0.101  0.842 


voltage  /  V  cr  voltage  /  V 


2720 


M.Jung,  K.A.  Williams  /  Journal  of  Power  Sources  196(2011)  2717-2724 


0  1  2  3  4  5 


Fig.  7.  Voltage  response  reproduced  by  Simulink®  using  the  calculated  parameters 
(a)  duty  ratio  is  50%  and  (b)  duty  ratio  is  15%. 


(open),  and  the  rate  is  high  when  the  switch  is  on  (closed),  so  the 
compressor  speed  generally  follows  the  input  signal.  However,  with 
increasing  frequencies,  the  compressor  speed  does  not  track  the 
input  signal  as  well.  This  is  believed  to  be  due  to  the  inertia  of  the 
compressor;  as  the  switch  is  closed  and  opened  more  frequently, 
the  compressor  cannot  speed  up  or  slow  down  quickly  enough  to 
track  the  varying  load.  At  the  frequency  range  between  about  1  Hz 
and  20  Hz,  the  period  of  the  compressor  operation  keeps  bounc¬ 
ing  back  and  forth  with  much  longer  frequencies  than  the  input 
frequency.  Eventually,  with  higher  input  frequencies,  the  compres¬ 
sor  speed  becomes  relatively  constant.  The  longer  frequency  of  the 
compressor  could  explain  the  discrepancy  between  experimental 
and  simulated  data  in  Fig.  8(a)  in  which  the  voltage  drop  of  the 
experiment  is  getting  bigger  due  to  the  continuous  operation  of 
the  compressor  of  which  the  magnitude  gradually  increases  for  a 
while.  The  results  at  higher  frequencies  in  Fig.  8(b)-(d)  show  dis¬ 
crepancies  during  early  stage  of  the  operation,  but  the  simulated 
results  matched  well  with  the  experimental  data  except  these  ini¬ 
tial  differences. 


3.2.  Analysis  of  dynamic  behavior 

The  experimental  data  can  show  only  the  total  voltage  drop,  but 
the  simulation  shown  in  Fig.  6  can  be  used  to  estimate  the  voltage 
drop  caused  by  each  RC  component.  The  simulation  results  of  each 
component  as  the  frequency  of  the  pulse  input  changes  are  plot¬ 
ted  in  Fig.  9.  In  the  result  for  the  charge  transfer  resistance  ( Rct\ 
the  full  voltage  oscillation  responding  to  the  input  signal  occurs  up 
to  about  20  Hz,  and  the  oscillation  amplitude  decreases  as  the  fre¬ 
quency  increases.  The  0.1  Hz  and  0.2  Hz  traces  are  not  visible  on  the 
plots  (Fig.  9(a)  and  (b))  because  they  are  underneath  the  1  Hz  trace. 
Due  to  the  small  value  of  Cdi  (0.01 -0.02  F)  and  Rct  (0.07-0.08  £2), 
there  are  still  oscillations,  even  at  the  relatively  high  frequency  of 
1 00  Hz.  Similar  results  were  obtained  for  other  two  RC  cells,  but  the 
frequency  ranges  at  which  oscillation  amplitudes  are  attenuated 
are  dependent  on  RC  values.  Although  different  values  for  the  down 
and  up  curves  were  used,  the  dynamic  behavior  of  each  RC  cells  are 
similar  between  both  curves.  The  results  are  summarized  in  Table  2. 
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Fig.  8.  Voltage  Merged  curve  of  down  and  up  curves  simulated  at  (a)  1  Hz,  (b)  10  Hz,  (c)  20  Hz,  and  (d)  100  Hz.  Modified  RLt  at  each  frequency  is  0.825,  0.805,  0.775,  0.775  £2 
for  the  down  curve,  and  0.842,  0.705,  0.692,  0.662  £2  for  the  up  curve  respectively. 


3.3.  Chemical  degradation  during  dynamic  operation 

As  described  previously,  chemical  degradation  continues  during 
non-operation  of  the  Ballard  Nexa  system.  Cessation  of  operation 
can  make  the  chemical  degradation  all  the  worse.  Therefore,  it  is 
important  to  ensure  how  long  the  fuel  cell  system  has  been  with¬ 
out  operation  before  performing  experiments  that  will  be  used 
to  assess  the  variation  of  membrane  degradation  as  a  function 
of  the  input  signal.  To  remove  the  effect  of  potential  chemical 
degradation  due  to  the  fuel  cell  system  not  being  used  for  peri¬ 
ods  of  time,  constant-current  aging  was  performed  on  the  Nexa 
system  before  each  such  experiment.  For  example,  since  the  fuel 
cell  system  had  not  been  operated  for  about  1  month  before  the 
first  experiment,  constant  current  aging  was  conducted  for  4.5  h 
for  the  first  experiment.  For  subsequent  experiments,  constant- 
current  operation  was  performed  for  2.5  h  to  mitigate  potential 
degradation  due  to  the  systems’  being  inactive  overnight.  More 
operating  details  are  provided  below  in  Table  3.  For  each  fre¬ 
quency,  dynamic  operation  of  the  fuel  cell  was  performed  for  a 
specific  duration.  At  fixed  times  during  operation,  effluent  sam¬ 
ples  were  taken  and  tested  for  the  concentration  of  HF.  For  all 
inputs  except  at  0.1  FIz,  test  durations  were  7  h.  For  the  0.1  Hz  input, 
an  additional  4.5  h  of  testing  were  performed,  to  provide  insight 
into  any  time-dependence  of  the  fluoride  emissions  over  that  long 
time  period.  No  such  dependences  were  found.  For  all  other  tests, 
the  duty  ratio  was  set  at  15%.  These  details  are  more  succinctly 
listed  below  in  Table  3,  and  the  results  are  shown  in  Fig.  10(a)  and 
(b). 


The  highest  level  of  fluoride  ion  released  in  effluent  in  the  first 
result  can  be  found  at  the  first  point  of  the  constant  aging,  suggest¬ 
ing  that  the  chemical  degradation  during  non-operation  influences 
the  fluoride  release.  From  the  fact  that  the  lowest  concentration 
of  fluoride  ion  in  the  Fig.  10(b)  is  at  the  first  point  of  the  dynamic 
operation,  it  is  thought  to  be  that  the  degradation  occurred  dur¬ 
ing  overnight  cessation  can  be  neglected  by  the  constant  aging  for 
2.5  h.  And  also,  the  fluoride  emission  increased  and  saturated  over 
about  8  h  during  the  dynamic  operation. 

The  frequency  effect  on  chemical  degradation  is  plotted  in 
Fig.  11.  Fluoride  emission  increases  up  to  peak  around  0.1  Hz  and 
then  decreases  with  further  frequency  increases.  In  examining  the 
data  in  Fig.  11(a),  it  appears  that  the  fluoride  emission  is  rela¬ 
tively  constant  for  frequencies  above  20  Hz.  If  it  is  assumed  that  the 
emissions  at  1  kHz  are  essentially  the  same  as  for  constant-current 
operation,  then  the  difference  in  fluoride  emissions  between  the 
constant  and  dynamic  operations  would  appear  to  be  on  the  order 
of  two,  rather  than  the  variation  of  a  factor  of  10  reported  by  Liu 
and  Case  [8].  This  could  be  from  different  experimental  conditions 
in  which  the  cyclic  operation  condition  included  open  circuit  [8]. 

It  is  interesting  to  compare  the  frequency  dependence  of  the 
fluoride  emission  with  the  voltage  response  of  the  first  RC  cell  (fti 
and  Ci )  listed  in  Table  2.  The  two  sets  of  data  are  both  shown  in 
Fig.  11(b).  In  that  figure,  the  frequency  range  can  be  divided  by 
three  zones:  a  full  oscillation  region  (DC-0.2  Hz),  a  transition  region 
(0.2-10  Hz),  and  no  oscillation  region  (>10Hz).  The  fluoride  emis¬ 
sions  increase  with  frequency  across  the  full  oscillation  region,  then 
decrease  across  the  transition  region  and  finally  become  indepen- 


Table  2 

Frequency  ranges  of  the  input  signal  leading  to  the  voltage  oscillation. 


Ret 

Ri  and  Ci 

R2  and  C2 

Full  oscillation  range 

DC-20  Hz 

DC-0.2  Hz 

DC-2  Hz 

Decreasing  amplitude  range 

>100  Hz 

0.2-10  Hz 

10  Hz-over  100  Hz 

No  oscillation  range 

- 

10  Hz-higher  frequency 

- 
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Fig.  9.  Voltage  response  dropped  by  Rct  simulated  using  the  parameters  from  (a)  the  down  curve  and  (b)  the  up  curve,  (c)  and  (d)  are  from  the  down  and  the  up  curves  by 
the  first  RC  cell  (Ri  and  C\ ),  and  (e)  and  (f)  by  the  second  RC  cell  ( R2  and  C2). 


Table  3 

Experimental  order  and  condition  for  dynamic  operation  of  the  Nexa  system. 


Day 

Freq.  of  input  signal 

Duty  ratio 

Cessation  period 

Constant  Current  aging 

Time  for  sampling  the  effluent 

1 

1  kHz 

15% 

~1  month 

4.5  h 

2  h,  4h,  5.5  h,  7h 

2 

0.1  Hz 

15% 

Overnight 

2.5  h 

2 h,  4h,  5.5 h,  7h,  8.5  h,  10h,  11.5h 

3 

0.0005  Hz 

15% 

Overnight 

2.5  h 

2  h,  4h,  5.5  h,  7h 

4 

1  Hz 

15% 

Overnight 

2.5  h 

2  h,  4h,  5.5  h,  7h 

5 

20  Hz 

15% 

Overnight 

2.5  h 

2  h,  4h,  5.5  h,  7h 

6 

0.001  Hz 

15% 

Overnight 

2.5  h 

2  h,  4h,  5.5  h,  7h 

7 

5  Hz 

15% 

Overnight 

2.5  h 

2  h,  4h,  5.5  h,  7h 

8 

0.01  Hz 

15% 

Overnight 

2.5  h 

2  h,  4h,  5.5  h,  7h 

dent  of  frequency  in  the  no  oscillation  region.  As  mentioned  in 
the  previous  discussion  of  modeling,  the  first  RC  cell  represents 
the  time  delay  corresponding  to  water  transport  in  the  membrane 
[11,12].  It  is  thus  supposed  that  the  accelerated  degradation  of  the 
membrane  (as  indicated  by  the  increased  fluoride  emissions)  is 
closely  related  to  water  transport  inside  the  membrane. 

4.  Discussion 

The  result  of  Fig.  1 1  showed  that  the  periodic  water  transport 
inside  the  membrane  plays  an  important  role  in  the  chemical  degra¬ 


dation.  In  general,  water  transport  in  the  membrane  is  considered 
as  to  result  in  dehydration  and  hydration  at  the  anode  and  cathode 
sides  of  the  membrane  respectively  because  of  electro-osmotic 
drag  [13].  There  are  four  kinds  of  water  movement  occurring  inside 
the  membrane,  electro-osmotic  drag,  back  diffusion,  evaporation 
into  air  at  the  cathode,  and  supply  from  external  humidifier  at 
the  cathode  and/or  the  anode  [13].  In  case  of  the  Ballard  Nexa 
system  having  very  thin  membrane  (3.5  x  10-5  m)  and  just  one 
side  humidifier  at  the  cathode  [14],  the  back  diffusion  of  the  water 
may  play  a  serious  role  in  the  water  movement  [13,15].  The  anode 
side  of  the  membrane  usually  experiences  the  humidity  cycle,  but 
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Fig.  10.  Fluoride  release  during  the  operation  (a)  at  1  kHz  and  (b)  at  0.1  Hz. 


the  cathode  side  is  fully  saturated  with  water  at  all  times.  There 
are  reports  indicating  that  humidity  cycles  may  induce  mechanical 
stresses  [16,17].  Moreover,  it  was  shown  that  membrane  thinning 
and  pin-holes  generation  were  accelerated  in  the  area  correspond¬ 
ing  to  the  stress  concentration  points  [18].  Therefore,  one  possible 
reason  for  accelerated  degradation  of  dynamic  operation  is  the 
periodic  stress  induced  by  the  humidity  cycles  even  though  the 
stress  is  not  enough  to  cause  mechanical  deformation. 

Another  possible  reason  suggested  in  this  paper  is  related  to 
diffusion  of  H2O2.  Electro-osmotic  drag  of  water  is  the  mass  flux 
caused  by  a  polarity  of  the  water  molecules  attracted  to  the  posi¬ 
tively  charged  protons  moving  from  the  anode  to  the  cathode  [19]. 
Since  H202  has  dipole  moment  of  2.2  Debye  which  is  bigger  than 
water  (1.85  Debye),  it  can  be  also  dragged  by  the  moving  protons 
and  periodically  transported  by  the  cyclic  operation  in  the  same 
fashion  as  the  water.  H2O2  is  generated  at  the  anodic  catalyst  layer 
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Fig.  11.  (a)  Fluoride  release  during  dynamic  operation  in  terms  of  frequency  vari¬ 
ation  and  (b)  relationship  between  the  voltage  oscillation  of  the  first  RC  cell  and 
fluoride  emission  in  terms  of  the  frequency  change. 


Distance 


Fig.  12.  Concentration  of  H202  inside  the  membrane  during  operation  (a)  at  open 
circuit,  (b)  during  switch  on  (large  current),  and  (c)  during  switch  off  (small  current). 


by  the  direct  combustion  of  hydrogen  with  crossover  oxygen  or 
incomplete  reduction  of  oxygen  on  the  Pt  surface  [1,6].  Further¬ 
more,  contaminants  such  as  chloride  ions  and  carbon  monoxide 
can  increase  the  yield  of  H202  drastically  [1].  Such  H202  forma¬ 
tion  reactions  can  also  occur  at  the  cathode  catalyst  layer  [1,20]. 
However,  H202  is  so  reactive  to  be  decomposed  into  water  rela¬ 
tively  quickly  on  the  Pt  surface  [20],  and  it  can  be  assumed  that 
there  is  a  little  chance  for  H202  to  be  accumulated  on  the  Pt  sur¬ 
face  before  it  diffuses  into  the  membrane.  This  assumption  can 
be  applied  to  just  the  cathode  catalyst  layer  because  a  monolayer 
of  hydrogen  on  the  Pt  surface  at  the  anode  inhibits  the  decom¬ 
position  of  H202.  For  this  reason,  it  is  widely  believed  that  H202 
formation  occurs  mainly  at  the  anodic  catalyst  layer  [  1 ,21  ].  With  no 
current,  H2O2  formed  at  the  anode  would  diffuse  into  the  cathode 
and  be  removed  by  decomposition  on  the  Pt  surface  or  evaporation 
into  the  air.  The  concentration  of  H202  at  the  anode  and  the  cath¬ 
ode  sides  of  the  membrane  would  depend  on  the  yield  of  H202  at 
the  anode  and  the  removal  rate  of  H202  at  the  cathode,  which  is 
described  in  Fig.  12(a).  Even  small  currents  could  reduce  the  H2O2 
concentration  at  the  anode  side  of  the  membrane  due  to  the  electro- 
osmotic  drag,  and  the  total  amount  of  H202  inside  the  membrane 
would  decrease  because  the  accumulation  of  H202  at  the  cathode 
could  expedite  the  removal  rate  as  shown  in  Fig.  12(b)  (a  bigger 
arrow),  which  can  explain  the  accelerated  degradation  at  open  cir¬ 
cuit  condition.  As  the  current  stops  or  decreases  during  the  dynamic 
operation  condition,  redistribution  of  H2O2  may  occur  inside  the 
membrane  due  to  the  back  diffusion  under  the  circumstance  of 
accumulated  H202  at  the  cathode  side  and  H202  generation  at  the 
anode  layer,  which  could  explain  deteriorated  degradation  during 
dynamic  operation.  In  addition,  the  serious  effect  of  the  low  humid¬ 
ification  on  the  degradation  mentioned  in  the  introduction  can  be 
explained  by  the  higher  dipole  moment  of  H202  than  that  of  the 


2724 


M.Jung,  K.A.  Williams  /  Journal  of  Power  Sources  196(2011)  2717-2724 


water.  It  was  reported  through  simulation  research  [22]  that  most 
water  molecules  are  constrained  by  sulfonate  heads  in  the  mem¬ 
brane  at  low  water  contents.  H202  would  also  be  constrained  by 
them,  besides,  more  strongly  than  water  due  to  its  larger  polarity, 
so  most  H202  could  make  contributions  in  formation  of  reactive 
oxygen  radicals  with  little  chance  to  diffuse  out  to  the  cathode. 

5.  Conclusions 

The  frequency  effect  of  the  dynamic  operation  on  the  chemical 
degradation  was  investigated  in  this  research.  Firstly,  the  Ballard 
Nexa  system  was  analyzed  by  a  modified  Randles  model  and  curve 
fitting,  and  dynamic  behavior  of  it  as  the  change  of  frequency  was 
simulated  using  MATLAB  Simulink®.  Secondly,  the  fluoride  release 
in  terms  of  the  frequency  change  was  evaluated  by  measuring  the 
concentration  of  fluoride  ion  in  the  drain  water  of  the  Nexa  system. 
The  comparison  between  above  two  results  showed  that  accelerat¬ 
ing  factor  of  the  dynamic  operation  on  the  chemical  degradation  is 
closely  related  to  water  transport  inside  the  membrane.  Two  pos¬ 
sibilities  of  the  accelerating  factor  were  presented.  One  is  the  stress 
induced  by  cyclic  humidity,  as  suggested  by  other  sources  and  the 
other  is  periodic  transport  of  H202  inside  the  membrane. 

Appendix  A.  Supplementary  data 

Supplementary  data  associated  with  this  article  can  be  found,  in 
the  online  version,  at  doi:10.1016/j.jpowsour.2010.10.105. 
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